Organization and instabilities of entangled active polar filaments. 
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We study the dynamics of an entangled, isotropic solution of polar filaments coupled by molecular 
motors which generate relative motion of the filaments in two and three dimensions. We investigate 
the stability of the homogeneous state for constant motor concentration taking into account excluded 
volume and entanglement. At low filament density the system develops a density instability, while 
at high filament density entanglement effects drive the instability of orientational fluctuations. 
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Cellular biology provides many realizations of pattern 
formation in dissipative nonequilibrium systems. An im- 
portant example is the collective behavior of the proteins 
that compose the cytoskeleton of eukaryotic cells. The 
cytoskeleton provides both the supporting structure of 
the cell and the vehicle for internal transport processes 
jij. It is a network of long protein filaments, mainly mi- 
crotubules, actin filaments and intermediate filaments, 
coupled by smaller proteins, such as molecular motors 
and cross-linkers. Motor proteins convert chemical en- 
ergy derived from the hydrolysis of ATP (Adenosine 
Triphosphate) into mechanical work, generating forces 
and motion of the filaments relative to each other in these 
active gels. 

Numerous in vitro experiments j^, |3|, ^, ^ have shown 
that mixtures of filaments and their associated mo- 
tor proteins self-organize into macroscopic symmetry- 
breaking structures, including radial arrays or asters and 
one-dimensional bundles. The nonequilibrium forces that 
give rise to these structures include the action of molec- 
ular motors and the polymcrization/depolymerization 
process of the filaments 0. Here we focus on the 
role of motor proteins and assume that the filaments 
have fixed length — a situation that can be achieved 
in vitro Q. A few analytical and numerical studies 
have investigated the emergence of these complex pat- 
terns g, H H |[ [fi], [□], [12). Continuum models of fila- 
ment/motor systems in two dimensions have been used 
to show that spatial patterns are obtained as nonequi- 
librium solutions of the system dynamics ||, [TT||. These 
models have ignored either filament diffusion or the 
motor action on orientational dynamics 0. A more mi- 
croscopic approach was taken by Kruse et al who con- 
sidered a dynamical model for the development of con- 
tractile and motile structures in one dimensional polar 
filament bundles, while igno ring steric and other interac- 
tions between the filaments Jh| [l^] . 

Many open questions remain concerning the role of the 
physical properties of the filament/motor gel in control- 
ling the formation of self-organized structures. Experi- 
ments have indicated that motor properties, such as their 



processivity - the fraction of time in a cycle a motor re- 
mains attached to the filament, strongly influence pattern 
formation. This is evident by comparing in vitro experi- 
ments in microtubules-kinesin to those in actin-myosin 
mixtures. At high motor concentration, microtubulc- 
kinesin mixtures readily organize in a variety of spatial 
patterns, provided the motors stall at the polymer ends 
before detaching g, In contrast, the homogeneous 
state is much more robust in the weakly coupled actin- 
myosin II systems, where spatially inhomogeneous struc- 
tures develop only upon depletion of ATP or at much 
higher filament concentration 13 . The physical charac- 



teristics of the filaments, such as their persistence length, 
may also contribute to the different behavior of these two 
active gels, as actin-myosin networks are more strongly 
entangled than microtubule- kinesin mixtures. Both an- 
alytic theories 10, [n| |l2| and simulations j|, || have so 
far entirely neglected entanglement. 

In this letter we start from a phcnomcnological model 
in the spirit of Kruse et al. |l2) and obtain a set of 
continuum equations to describe the dynamics and or- 
ganization of polar filaments driven by molecular motors 
in an unconfined geometry in (quasi-)two and three di- 
mensions (d = 2,3). By modeling the motor-filament 
interaction microscopically, we can calculate the magni- 
tude and, most importantly, the sign of the parameters 
of the continuum equations, which cannot be obtained by 
symmetry arguments. We consider a isotropic filament 
solution and include the effects of entanglement and ex- 
cluded volume. Our result is a phase diagram (Fig. |^) as 
a function of the filament density and motor properties. 

The filaments are modeled as rigid rods of length I 
and diameter b « I. Each filament is identified by the 
position r of its center of mass and a unit vector n point- 
ing towards the polar end. Taking into account filament 
transport, the normalized filament probability distribu- 
tion function, ^(r,n,t), obeys a conservation law [ pd| , 

a t * + V-J + ^-J r = 0, (1) 

where 1Z = n x d& is the rotation operator. The transla- 
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tional and rotational currents J and J r are given by 

Ji = -Diidi* - * d 3 V cx + Jf ct , (2) 

Jl = -D r Ki* - p- * n t v^ + j; /act , (3) 

where i = 1, and Dy = Duhifij + D±(8ij — hifij) is 
the translational diffusion tensor and D r the rotational 
diffusion constant. The potential V cx incorporates ex- 
cluded volume effects that play an important role in sta- 
bilizing time-dependent solutions. It is given by fcgT 
times the probability of finding another rod in the inter- 
action area of a given rod, 

V cx (r, fix) = k B T [ [ (r + £, n 2 ) , (4) 

where the prime restricts the integral to the interaction 
volume, corresponding to the region where the two fila- 
ments touch at at least one point. The volume of this 
region is Vint = Vo^/l C&1 • n 2 ) 2 , with u = l 2 b d ~ 2 and 
l 2 \J\ — (ni • ii 2 ) 2 > b 2 . The active currents are given by 

J act (r, *i)=[ / v(C, n 1; n 2 )*(r, n 1 )*(r + {, n 2 ), (5) 

J r / act (r, n x ) =/ fix, fi 2 )¥(r, n x )*(r + £, n 2 ), (6) 

where v = £ and w = ni — n 2 are the relative linear 
and angular velocities of two filaments, with the dot de- 
noting a time derivative. The model naturally contains 
two competing dynamics. The first is the diffusion of 
hard rods, which at high density must include excluded 
volume and entanglement. The second is the local driv- 
ing force coming from the interaction with the motors. 
This depends on the polarity of the filaments and breaks 
the n — > — n symmetry of the hard rod fluid, allowing for 
states of broken symmetry, where the filaments acquire 
a nonvanishing mean orientation. 

In the absence of external forces and torques the to- 
tal linear and angular velocity of an interacting pair are 
conserved. This requires v(£, ni,n 2 ) = — v(— £, fii, n 2 ) 
and u>(£, ni,n 2 ) = — u:{— £, &i, n 2 ). Rotational 
and translational invariance requires v(£, ni,n 2 ) = 
-v(-£, -fix, -n 2 ) and u>(£, rii, n 2 ) = <*>(-£, -h 1 , -n 2 ). 
The simplest form of the linear and angular velocities can 
be written as 

v = a_ g(j_+ni -n 2 ) [3 n 2 - n x 

V ~ 2Z ^1 - (fii • n 2 ) 2 2 ^1 - • fi 2 )2 ' 



v/i-(fii-ft 2 ) 2 

The velocities have been normalized with the volume of 
interaction. The parameters a, [3 and 7 are the rates 
for the various motor-induced translations and rotations 
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FIG. 1: Cartoons of various motor-induced filament inter- 
actions. All interactions are viewed from the rest frame of 
filament 2. The initial and final position of filament f are 
shown as a thick and a dashed arrow, respectively. The zig- 
zaged lines represent motors, (a) The contribution to v pro- 
portional to a which is along the direction of the relative dis- 
placement £ of the centers of mass of the two filaments (thin 
arrow). Motors also drive rotation of filament f as indicated. 
The contribution to v proportional to ji is illustrated in (b) 
and (c) for two filaments with £ = and ni • n2 > (b) and 
ni ■ 112 < (c). In both cases the translation at a rate ft in 
the direction of n2 — ni (thin arrow) tends to bring the polar 
heads of the two filaments to the same spatial location. In 
(b) the counterclockwise rotation aligns the filaments, while 
in (c) the clockwise rotation anti-aligns and separates. 

(Fig. |l|). The contribution proportional to a depends on 
the separation of the centers of the filaments and results 
from a difference in motor activity between the ends and 
mid-points of the filaments. It tends to align the cen- 
ters of mass and polar heads of the filament pair (see 
Fig 2(a)). The contribution proportional to (3 vanishes 
for aligned filaments and can separate antiparallel fila- 
ments, as illustrated in Fig. 2(c). This mechanism yields 
both translational and rotational currents. The prefactor 
(rii ■ n 2 ) in the angular velocity guarantees that motors 
preferentially bind to two filaments that are at an angle 
smaller than 7r/2. The 7 term has no effect on perpen- 
dicular filaments. Assuming uniform motor density p m , 
from simple mechanical models of motors M, we esti- 
mate a ~ (3 ~ 7Z ~ pmlb 4>(s c /t c ), with s c the motor 
step length per cycle, r c the time for one cycle and <p the 
duty ratio. 

We are interested in describing the dynamics of ac- 
tive filaments on length scales large compared to the fil- 
aments size, /. We can then expand the concentration of 
filaments $>(r + £, n 2 ) near its value at r, 

*(r + £, n 2 ) = *(r, n 2 ) + £„e„ ■ W(r, n 2 ) 

+ l^ m (e n ■ V)(e m • V)*(r, n 2 ) + 0(f). (9) 

We have introduced a set of orthogonal unit vectors, 
(ei,e 2 ,z), that provides a natural coordinate system for 
the problem. The unit vector z is normal to the plane 
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passing through the point of contact of the two filaments 
and containing the unit vectors fii and ii 2 . The vectors 
ex = (fix + n 2 )/|nx + n 2 and e 2 = signfnx • n 2 )(n 2 - 
ni)/|n 2 — ftxl are orthogonal unit vectors in this plane. 
Neglecting the out-of-plane separation (of order b) be- 
tween the centers of mass of the two filaments, the vector 
£ is written in this coordinate system as £ = £ rl e„, where 
summation over n = 1,2 is intended. We assume that 
on large scales the filament dynamics can be described in 
terms of the filaments density p(r) and the local filament 
orientation t(r) defined as the first two moments of the 
distribution \P(r, n, t), 



t(r, t) 



du ( „ )*(r, n,i) 



(10) 



Coarse-grained equations for p and t can be obtained by 
inserting Eq. (|^) in the expressions for the active cur- 
rents and for V ex , writing the density \&(r, n, t) in the 
form of an exact moment expansion, and retaining only 
the first two moments in this expansion. For brevity, 
we only display here the dynamical equations linearized 
about a homogeneous state, with constant density po and 
an isotropic orientational distribution of filaments, cor- 
responding to t = 0. The full and rather cumbersome 
nonlinear equations will be given elsewhere 15 . Letting 
p = po + Sp and keeping only terms up to third order in 
the gradients, the linearized equations are given by 



d t 5p = ~ [D l{ + (d - l)Dx] (1 + v oP o)V 2 5p - 



advoPo 
12d 



V 2 Sp 



f3l 2 v oPo (2d+l) 
24d(d + 2) 



V 2 (V-t) , 
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dtU 



1 



' ' d + 2 
alvopo 



[(d+l)L> ± +D||]V 2 ^ 



12d(d + 2) 



[V 2 U + 2ftV • t] 



d + 2 
'-di5p 



v oPo as , I3l 2 v p (2d+1) 



£>±)ftV-t 

d l V 2 dp . 



24d 2 (d + 2) 



(12) 



The local orientation is not a conserved variable and de- 
cays at a rate ~ D r . Both equations display the com- 
petition of diffusive terms (oc -DV 2 ) and pattern- forming 
terms (oc — aV 2 ). The linear instability of the homoge- 
neous state occurs when the pattern-forming terms dom- 
inate. To linear order, the contribution from the rota- 
tional current vanishes and excluded volume corrections 
only appear in the density equation. 

We now study the stability of the homogeneous state. 
We expand the fields in Fourier components, 5p(r) = 
Ek/°k eJk r an d t(r) = X)k t ke lk r , and separate tk in 
its component longitudinal and transverse to k, namely 
tk = k • tk and t£ — k x tk, with k = k/|k|. In d 
dimensions there are d — 1 degenerate transverse modes 
describing the decay of fluctuations in t£, with rate 



A T (*0 



D r 



d+2 



(d+l)D A 



■Du 



alvoPo 
12d 



(13) 



There are two coupled modes describing the decay of den- 
sity and i£ fluctuations, given by 

A±(*) = l{ M n + M 22 ± V(Mn - M22) 2 + AM 12 M 2 i } , 

(14) 

with 



M n =~ +(d- l)D ± )(l + v QPo ) alvoPo 
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M- 



22 



■D r 



d + 2 



3Du + (d-l)D A 



alv po 



Ad 



M 



ik 



3 /3l 2 v p 2d+l 



24 d(d + 2) ' 
., ., JvoPo l 2 k 2 2d+l 

M 2 l = tk : 1 



d \ 24 d(d + 2) 



(15) 



At long wavelength, A + = A p vanishes as k 2 and describes 
the decay of density fluctuations, while A_ = is a 
kinetic mode (i.e., it is finite as k — > 0) and describes 
the decay of longitudinal orientation fluctuations. The 
hydrodynamic density mode goes unstable on all length 
scales for a > a { P ~ (1 + (d - 1)D X /Z>||)[1 + (/5 ) -1 ], 
where a = al/(12Du) and p = PqVq. 

Dilute Solutions. For dilute solutions of long thin 
rods the diffusion constants are D± = Dn/2 = D/2 and 
D r = 6D/1 2 , with D = k B T1n(l/b)/(2irr)l) and 77 the sol- 
vent viscosity H]. In this regime the instability of the 
homogeneous solution occurs at and is associated 
with the density mode. 

Semi-dilute Solutions. For semidilute solutions, the 
dynamics is modified by the topological constraint that 
the filaments cannot pass through each other. This 
constraint can be modeled by a tube lead- 
ing to modified transverse and rotational diffusivities, 
Dj_ ~ D/[2(l + cj_p {l/b) d - 2 ) 2 } and D r ~ 6L»/[^ 2 (1 + 
CrPo(l/b) d ~ 2 ) 2 ) with c r t ± numbers of order unity and Dm 
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FIG. 2: Linear modes and phase diagram {a, po} showing the 
phase boundaries aa 1 , Qc between the homogeneous state 
and the density and orientationally inhomogeneous state. The 
modes on the right corresponds to al/D = 10, po = 10, l/b = 
100 while those on the left al/D = 1900, p = 1, l/b = 100. 



essentially unaffected by entanglement |14|, M. The ki- 
netic modes A_ and At describing the decay of orien- 
tation fluctuations can then become unstable before the 
density mode. This instability occurs at a c ~ 1 / po for 

k > ko, with fco — (X/lpo)\l w^jizz^- A stability diagram 



\/l 3 /£ p leading to vq 



/"HUjy 3(&p -d)- 

in the (a, po) plane is shown in Fig. |2|. 

Semiflexibility. Semiflexible filaments of persistence 
length £ p > I can be modeled as fuzzy rods with 
a diameter given by b 

j2+3(d-2)/2£-(<2-2)/2 

Discussion. The linear instability of the homoge- 
neous state is controlled by the parameter a that drives 
filament bunching, while f3 and 7 play no role |17| , At 
low filament densities, well below the critical density for 
the nematic transition, density fluctuations become un- 
stable on all length scales at ctp , signaling the onset of 
a state with inhomogeneous density. For p > p c ~ 2, 



where the critical curves and olc' cross, the instabil- 
ity occurs at the lower value a c and it corresponds to 
a short-scale instability of orientation fluctuations, while 
the density remains homogeneous. In this regime the fil- 
ament solution has a substantial degree of short range 
nematic order and both transverse and rotational diffu- 
sion are strongly impeded. The orientational degree of 
freedom, t, becomes essentially conserved as D r is van- 
ishingly small. The mechanism for the orientational in- 
stability of dense filament solutions obtained here is dis- 
tinct from that proposed by Lee and Kardar [jll| . These 
authors modeled the dynamic in the nematic phase and 
incorporated an inhomogeneous motor density. In their 
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model the instability is controlled by the nonlinear cou- 
pling of orientation to motor density fluctuations. Here, 
in contrast the instability occurs to linear order and is 
driven by the pattern- forming terms in Eq. (p^). 

A number of open questions remain. It will be rela- 
tively straightforward to include the motor transport us- 
ing our formalism. This is important for very processive 
motors and at low motor densities. An inhomogeneous 
motor density may also be required for the formation of 
stable asters and vortices at low filament concentration. 
It will be also very interesting to ask about the response 
of the motor/filament system to shear. This work will be 
presented elsewhere Jl8| . 
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